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Abstract 

We demonstrate the use of an ordinary white-light microscope for the study of the q-dependent 
dynamics of colloidal dispersions. Time series of digital video images are acquired in bright held 
with a fast camera and image differences are Fourier-analyzed as a function of the time delay 
between them. This allows for the characterization of the particle dynamics independent on whether 
they can be resolved individually or not. The characteristic times are measured in a wide range 
of wavevectors and the results are found to be in good agreement with the theoretically expected 
values for Brownian motion in a viscous medium. 
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Microscopy and light scattering are widely nsed in physics, chemistry, biology and medical 
laboratories to access information on the strnctnre and dynamics of mesoscopic systems. 
While microscopy gives direct access to real space images, scattering techniqnes work in 
reciprocal space, where information on the strnctnre and dynamics of the system is obtained 


Q. 


respectively from the angnlar and time dependence of the scattered light intensity Hi. These 
two complementary techniqnes have in general very different experimental reqnirements. 
White light sonrces are nsnal choices in microscopy, while a certain degree of coherence of 
the illnminating beam is reqnired in scattering experiments; this is usnally achieved by nsing 
a laser. In the past many attempts have been made to bnild a scattering apparatns based on a 


microscope; all o 
example Refs. 


these attempts involved the use of a laser as an illumination source (see for 
6|] and references therein). In some cases special care was taken to ensure 
the capability to perform microscopy and scattering experiments simultaneously, thereby 
allowing for a powerful combination of the complementary information obtained by both 


techniques 




More recently, microscope-based laser Dynamic Light Scattering (DLS) 


experiments have been developed to study the dynamic properties of samples of biological 


interest, such as living macrophage and red blood cells {^,0]. In practice, due to the intrinsic 
difficulties in building such instruments, the use of such techniques has been restricted to 
those laboratories, where a sufficient expertise in the realization of optical instrumentation 
was at hand. 

In this letter, we present a conceptual scheme to interpret and analyze microscopy images 
that are obtained from samples containing moving entities. This technique, which we term 
Differential Dynamic Microscopy (DDM), does not entail any special experimental require¬ 
ments, being based on the use of a standard light microscope with a normal illumination 


source and a digital video camera. By using the tools of Fourier Optics we provide 
the means to access information about the sample dynamics that are equivalent to the one 
obtained in multi-angle dynamic light scattering (DLS) experiments [Sj. We test DDM 
by analyzing time sequences of microscopy images obtained with an aqueous dispersion of 
colloidal particles with diameter 73 nm, well below the resolution limit of the microscope. 
Additional tests are performed with larger particles (diameter 420 nm). Our results are 
in good agreement with the predictions from the theory of Brownian motion in a viscous 
medium, thereby validating the proposed analysis. We believe that DDM is an interesting 
complement to the well established DLS 9| and microscope video tracking techniques. 
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Our experimental setup simply consists in an inverted microscope (LEICA DM IRB) 
equipped with a complementary metal-oxide-semiconductor (CMOS) camera (IDT X-Stream 
XS-3, 1280 X 1024 pixels, pixel size=12/im). No modifications are made to the microscope. 
The sample is illuminated by focusing white light on the sample with a condenser lens 
(numerical aperture 0.9); for detection we use a standard objective (magnification M = 63x, 
numerical aperture NA = 0.7). The effective pixel size, accounting for the magnification M, 
is then given by dpix = 0.19/rm, value smaller than the microscope resolution limit imposed 
by diffraction. A capillary tube with rectangular section (Vitrocom, Inc.) is filled with an 
aqueous dispersion of polystyrene spheres of diameter 73 nm (Duke Scientific, part number 
3070A). The concentration is 1% by weight fraction, sufficiently low to neglect the effect 
of interactions between individual spheres. The thickness of the capillary tube along the 
microscope optical axis is lOO/rm and a region at the center of the capillary tube is imaged 
onto the camera sensor. 

A sequence of A^ = 1000 images is acquired with a sampling rate of 100 samples/s and 
with an exposure time of 1 ms. In Fig.1(a) we show a typical image: the image appears 
as a bright background with some dark spots. These dark spots are dust particles located 
mainly on the protective glass of the camera sensor. The signal due to the particles is barely 
visible as its contribution is small compared to the large background signal. However, the 
background is a static contribution, while the Brownian motion of the particles causes an 
incessant rearrangement of their configuration in space. By subtracting two images grabbed 
at different times we can eliminate the static contributions in our image and bring the 
displaced particles to light. Examples of difference images obtained by this subtraction 
procedure are shown in Fig.l(b) and 1(c); (b) and (c) are respectively the difference images 
obtained by subtracting two images separated in time by At = 10ms and At = 200ms. 
Clearly, the grainy appearance of the difference images increases with increasing time delay 
At, revealing the increasing displacement of the particles. 

We can quantify this increase by calculating the spatial variance 

(T^(At) = J \D{x,y; At)f dxdy, (1) 

of the difference signal D{x,y, At) = I{x,y,At) — l{x,y;0) as a function of At. Here, 
I{x,y,t) is the intensity obtained in the sensor plane {x,y) at the time t, where we define 
the acquisition time of the first image to be t = 0 s. The actual choice of the reference 
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FIG. 1: (a) Microscope image of an aqueous dispersion of 73 nm particles at a weight fraction of 
1%. The particles are below the resolution of the microscope and the signal generated by them is 
very small, (b) Result of the subtraction of two images taken 10 ms apart in time. The signal due 
to the particles is now visible, (c) Result of the subtraction of two images 200 ms apart in time. 
The signal has increased. The size of each panel corresponds to 47 fim in the sample. 


image is not crucial in the case of Brownian motion, since the statistical properties of the 
investigated dynamics do not change in time. This can be appreciated in the movie clips 
shown in [h|, where the time delay between two subtracted images is kept constant at 
At = 200 ms. Consequently, for every At we calculate the variance cr^(At) by averaging 
over many statistical realizations of the signal, typically 100, to improve the statistical 
accuracy. As shown in Fig.2, the variance a‘^{At), increases very rapidly with At and 
attains a steady state plateau at long times, when all the particles have moved far away 
from their initial position, such that the two images used to obtain the difference image 
are statistically uncorrelated. This dependence of the variance on At clearly indicates that 
difference images can be used to monitor the particle dynamics. The question then arises 
whether it is possible to extract quantitative information on the dynamics of the system 
from the difference images. The answer is yes, but its justihcation requires the use of two 
steps. 

The hrst one is to consider that the variance, dehned in EqU] can be calculated in 
the Fourier space as well. If we dehne the 2D Fourier transform of D{x,y]At) as 
^{ux, Uy] At) = / D{x, y] At)exp [—i27i {uxX + Uyy)] dxdy, the use of the Parseval theorem 
[TJ guarantees that 

a^{At) = J \FD{ux,Uy-, At)f duxduy. (2) 
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FIG. 2: The particles rearrangement due to their Brownian motion causes an increase of the 
variance cj^(At) of the difference signal D{x, y; At). With increasing time delay the two subtracted 
images become progressively uncorrelated, such that cr^(At) saturates at large At. 


The equivalence of Eq{T] and [2] shows that the total energy content of D in real space is the 
same as the one of in Fourier space. However, in Fourier space it is possible to isolate 
every Fourier component and study the growth of its amplitude as a function of At. 

To clarify whether this mode decomposition provides insight about the dynamics of the 
system we need to go through the second step. This step consists in realizing that a Fourier 
component of the microscope image is associated in a simple way to a scattering angle and in 
turn to a scattering wavevector. According to Abbe’s theory of microscope image formation 


3, 


12j, the object to be imaged on the camera sensor can be thought as the superposition 


of different Fourier components characterized by a spatial frequency Uobj- Fach frequency 
component of the object acts as a grating and diffracts light at an angle t} = sin~^{Xuohj) 
with respect to the microscope optical axis, where A is the wavelength of light. On the camera 
sensor, the diffracted plane wave will produce a sinusoidal fringe pattern, characterized by a 
spatial frequency Udet = sin['d)/\ = Uobj, independent of the wavelength A. This guarantees 
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FIG. 3: Growth of |FD(g;At)p with At for three values of q: q = 2.6/im“^(black squares) , 
q = 3.9/um“^(red up-triangles) and q = 4.5/im“^(blue down-triangles). The continuous lines are 
fits of the data to EqlSl 

that we can associate in a nniqne way, a spatial freqnency Udet on the detector to a scattering 
angle 'd. 

To investigate the angular dependent dynamics of our system we thus use the Fourier 
power spectrum \F£,{ux,Uy] At)f, which we simply calculate by applying a Fast Fourier 
Transform to the difference images. For isotropic samples, like the one investigated, the 
expectation value of \F]:){ux,Uy] At)f is rotationally invariant in the {ux,Uy) plane. This 
allows us to perform azimuthal averages and to treat the one-dimensional power spectrum 
|Fd(m; At)|^, where u = + Uy. For easy comparison with scattering experiments we 

use wavevectors q = 27iu instead of spatial frequencies u. Processing the data at several At 
enables us to report the one-dimensional power spectrum as a function of the 

time delay, as shown for three values oi q {q = 2.6/im“^, q = and q = 4.5/im“^) 

in Fig.3. As for the variance a‘^{At) we find that |Fn)(?; At)f increases rapidly with At to 
then saturate at longer time delays. The cross-over to the plateau value systematically shifts 
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to smaller At as the wavevector q increases, denoting a g-dependence in the characteristic 
time of the system. 

For a quantitative description of this behavior we consider that fluctuations in the inten¬ 
sity of the original images are due to concentration fluctuations in the sample. For Brownian 
diffusion it is well known that every Fourier concentration mode decays exponentially in time, 
exp(—At/r(g)), with a characteristic time r(g) = where Dm is the mass diffusion 

coefficient of the particles Iq]. For the analysis of the difference images it is then easy to 
show that 


iFniq; Af)|^ = A{q)[l - exp{-At/T{q))] + B{q). 


( 3 ) 


This equation contains two main contributions: the term B{q) is due to the power spec¬ 
trum of the camera noise and is present even in absence of the particles. The term 
A{q)[l — exp{—At/T{q))] describes the contribution associated to the particles. A quantita¬ 
tive description of the term A{q) requires the knowledge of the exact relationship between 
intensity and concentration fluctuations, which is beyond the scope of this paper. However 
once q is hxed, we can treat A{q) and B{q) as simple adjusting parameters and extract 
the characteristic time r(g) by simply considering the Af-dependence of |iA)(g; Af)p. This 
is very similar to what has been done in Ref. where an expression similar to Eqjn] 

has been proposed to analyze sequences of shadowgraph images for the characterization of 
nonequilibrium concentration fluctuations in a binary mixture. 

Thus, we £t our data to Eqj3] to obtain r(g), which we report in Fig.4 for q ranging 
from 0.4 to 5 pm~^ (black open circles). On the same graph the dotted black line 

denotes the expected behavior calculated by using the formula r(g) = IjDmq^i where Dm 
is estimated using the Stokes-Einstein relation Dm = ksT/(STrpd) = 5.98 pm?/s] ks is the 
Boltzmann constant, T the absolute temperature, p the solvent viscosity and d the diameter 
of the particles. The experimental results are in good agreement with the theoretically 
expected values over more than one decade in q, showing that DDM is a reliable means to 
determine the dynamical properties of the sample. To our knowledge this result provides the 
hrst experimental evidence that white light, bright-held microscopy can be used to monitor 
the g-dependent dynamics of particles. 

We also perform experiments by using particles with a diameter of 420 nm, value com¬ 
parable to the point spread function of the microscope. In Fig.5 we present two images of 
this sample separated in time by 10 s. In contrast to the previous case the particles are 
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FIG. 4: Characteristic decay time r versus the wavevector q for the 73 nm (black circles) and the 
420 nm (blue squares) particles. Fitting of the data gives respectively the values Dm = 6.2 ± 0.3 
[irn? js and Dm = 1.1 ± 0.2 [iw? js. The lines are theoretical predictions calculated by using the 
Stokes-Einstein relation with no adjustable parameters. 


fairly visible, though the static signal due to dust along the optical path still represents a 
significant contribution. However, the subtraction procedure described above can still be 
applied to visualize the particles motion and to eliminate the static noise, as shown in the 
movie clip presented in [ll|. The image sequence is analyzed as described above and the 
results are presented in Fig.4 as blue open squares. A good agreement with the theoreti¬ 
cally expected behavior [Dm = 1.04/im^/s, continuous blue line) is also found in this case, 
demonstrating that DDM is capable of monitoring the dynamics of particles that can be 
either resolved individually or not. This suggests that DDM could be a valid complement to 
tracking techniques for relatively big particles, with the advantage of being also appli¬ 
cable to concentrated systems, where tracking becomes difficult, and to fluctuating systems 
which are not composed by particles. Typical examples are membranes, interfaces or liquid 
crystals. For small particles DDM can be subsidiary to DLS experiments Qj, because it can 
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FIG. 5: Microscope images taken 10 s apart. The side of each panel corresponds to 47//m in the 
sample. The sample is a colloidal dispersion of 420 nm particles. The more visible objects are dust 
particles on the CMOS sensor which are also present in Fig. 1(a), grabbed with another sample. 
The particles are the less visible objects that change position between image (a) and (b). The 
smearing effect due to the finite point spread function of the microscope is visible. 


operate at small wavevectors where stray light makes the use of DLS notoriously difficult. 
The range of accessible time-scales is limited by the inverse of the acquisition frame rate on 
the lower side and by the overall duration of the experiment on the higher end. 

A possible extension of this work could be the investigation of dynamically heterogeneous 
systems, where a space resolved investigation of the dynamics is of fundamental importance 
15[. Our technique could also be applied to absorbing samples and to fluorescence mi¬ 
croscopy, which is widely used for the investigation of biological samples 16|. Finally, we 
believe that DDM can be prohtably used to perform microrheology experiments as pioneered 
in Ref. Q- 

The authors thank M. Carpineti, F. Giavazzi and A. Vailati for stimulating discussions 
and for critical reading of the manuscript. T. Bellini, P. Cicuta, F. Croccolo, M. Giglio, R. 
Piazza, G. Takacs and D. Weitz are thanked for fruitful discussions. R.G. acknowledges h- 
nancial support from the European Union (Marie Gurie Intra-European Fellowship, contract 
EIF - 038772) and V.T. hnancial support from the Swiss National Science Foundation. 


9 








* Electronic address: roberto.cerbino@unimi.it 

[1] B. Chu, Laser Light Scattering: Basic Principles and Practice, (Dover, New York, 2007, 2nd 
ed) 

[2] P. D. Kaplan, V. Trappe, and D. A. Weitz, Appl. Opt. 38, 4151 (1999). 

[3] M. T. Valentine, A. K. Popp, P. D. Kaplan, and D. A. Weitz, Opt. Lett. 26, 890 (2001). 

[4] R. Dzakpasu, and D. Axelrod, Biophys. J. 87, 1279 (2004). 

[5] R. Dzakpasu, and D. Axelrod, Biophys. J. 87, 1288 (2004). 

[6] M. S. Amin, Y. Park, N. Lue, R. R. Dasari,K. Badizadegan, M. S. Feld, and G. Popescu., 
Opt. Exp. 15 , 17001 (2007). 

[7] J. W. Goodman, Introduction to Fourier Optics, (Roberts&Gompany, Englewood, 2005, 3rd 
ed). 

[8] L. Gipelletti, and D. A. Weitz, Rev. Sci. Instrum. 70 , 3214 (1999). 

[9] B. J. Berne, R. Pecora, Dynamic Light Scattering: With Applications to Chemistry, Biology, 
and Physics, (Dover, New York, 2000). 

[10] J. C. Crocker, and D. G. Grier, J. Colloid Interface Sci. 179, 298 (1996). 

[11] See EPAPS Document No. [number will be inserted by publisher] for 50 Hz movie clips compar¬ 
ing the time evolution of the original images to the one of difference images with At = 200ms. 

73nm.mpg and 420nm.mpg are for 73 nm and 420 nm correspondingly. For more information 
on EPAPS, see http://www.aip.org/pubservs/epaps.htmL 

[12] H. Kohler, Opt. Acta. 28, 1691 (1981). 

[13] R. N. Bracewell, The Fourier Transform and Its Application, (McGraw-Hill, New York, 1999, 
3rd ed). 

[14] F. Croccolo, D. Brogioli, A. Vailati, M. Giglio and D.S. Cannell, Ann. N.Y. Acad. Sci. 1077 , 
365 (2006). 

[15] L. Gipelletti, and L. Ramos, J. Phys. : Condens. Matter 17 , R253 (2005). 

[16] J. W. Lichtman, and J. A. Conchello, Nat. Methods 2, 910 (2005). 

[17] T. G. Mason, and D. A. Weitz, Phys. Rev. Lett. 74, 1250 (1995). 


10 




